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ABSTRACT 



We observed VB 10 in August and September 2009 using the FORS2 camera of the VLT with the aim of measuring its astrometric 
motion and of probing the presence of the announced planet VB 10b. We used the published STEPS astrometric positions of VB 10 
over a time-span of 9 years, which allowed us to compare the expected motion of VB 10 due to parallax and proper motion with the 
observed motion and to compute precise deviations. The achieved single-epoch precisions of our observations are about 0.1 mas and 
the data showed no significant residual trend, while the presence of the planet should have induced an apparent proper motion larger 
than 10 mas yr -1 . Subtraction of the predicted orbital motion from the observed data produces a large trend in position residuals of 
VB 10. We estimated the probability that this trend is caused by random noise. Taking all the uncertainties into account and using 
Monte-Carlo resampling of the data, we are able to reject the existence of VB 10b with the announced mass of 6.4 Mj with the false 
alarm probability of only 5 ■ 10~ 4 . A 3.2 Mj planet is also rejected with a false alarm probability of 0.023. 

Key words, astrometry - technique: high angular resolution - planetary systems - stars: individual: VB 10 



1. Introduction 

The detection of a planetary-mass companion to the nearby 
M8 ultracool dwarf VB 10 (GJ 752B) was announced by 
Pravdo & Shaklan (2009). Using ground-based astrometric ob- 
servations with the STEPS camera of the 5 m Palomar telescope 
obtained over nine years, the authors derived a full astrometric 
solution of the system that suggested the presence of a planetary 
companion with a 271 day orbital period and a 6.4 Mj mass. 
Even though the VB 10b astrometric signature is large with a 
peak-to-peak value of 8 mas, the detection is debatable due to 
imperfect orbital sampling and comparably low astrometric ac- 
curacy of the individual measurements. 

Zapatero Osorio et al. (2009) presented high-resolution 
near-infrared observations of VB 10 obtained between 2001 and 
2008 with the NIRSPEC instrument on the Keck II telescope. 
They emphasized on the necessity to have a better sampling of 
the orbital phase to precisely constrain the orbital parameters 
and the individual masses of the system. Bean et al. (2010) used 
the near-infrared CRIRES spectrograph and did not detect the 
planet's signature in radial-velocity (RV). Anglada-Escude et al. 
(2010) performed a joint analysis of precise RV observations 
over 175 days with MIKE/Magellan and ESPaDOnS/CFHT 
spectrographs, RV measurements by Zapatero Osorio et al. 
(2009), and astrometric data by Pravdo & Shaklan (2009). They 
showed that the observed astrometric motion is not due to an 
unseen companion. 

We present results of the astrometric search for the planet 
using the FORS2 camera of the VLT, the only instrument which 



* Based on observations made with ESO telescopes at the La Silla 
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can achieve astrometric precision of individual measurements of 
0.05-0.1 mas (Lazorenko et al. 2009). Despite the short obser- 
vation period of two months, we obtained data sufficiently good 
to verify the hypothesis on the possible existence of VB 10b. 



2. Observations 

Observations were made with the FORS2/VLT camera 
(Appenzelleret al. 1998) in imaging mode during five nights 
from 2 August to 25 September 2010. The images cover a 
4.2 x 4.2'field of view with a pixel (px) scale of ~0.1"/px. For 
each night, we obtained 21 to 72 frames of 7 s exposures with 
the /^special filter. Seeing varied from 0.55"to 0.9". To reduce dif- 
ferential chromatic refraction (DCR) of the atmosphere, obser- 
vations were made near meridian. 

Because the observed motion of VB 10 in the sky is large 
(over 2 px), we had to know the precise pixel scale to convert 
pixel to arcsec units without loss of accuracy. From CMC 14 and 
NOMAD catalogue positions of stars in our reference frame, we 
derived acceptably precise scales of 0.12538 +0.00007 "/px and 
0.12567 + 0.00005 "/px in RA and Decl, respectively. 



3. Data reduction and analysis 

Our null hypothesis was that the planet VB 10b exists. To test 
this hypothesis, we computed the residuals of our measurements 
compared to the position of VB 10 expected from the combina- 
tion of proper, parallactic, and orbital motions. By determining 
the probability that the measured residuals are compatible with 
the expected mean value of zero, we could confirm or reject the 
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null hypothesis. We also performed the corresponding computa- 
tions assuming that the planet does not exist. 

In all cases, we considered that the astrometric signal *P 
in VB 10 position induced by the planet is defined by the 
orbital elements of VB 10b given in the discovery paper by 
Pravdo & Shaklan (2009). The predicted orbital motion over the 
measurement timespans of 17 and 54 days is 0.48 and 1.95 mas 
in RA and 0.92 and 3.28 mas in Decl, respectively. Compared 
to the astrometric precision of FORS2 of about 0.1 mas per 
epoch, these displacements are large and should be detected in 
this study. 

3.1. Photocentre determination 

Raw images were flat-fielded and bias-subtracted to exclude 
pixel-to-pixel variations of the CCD sensitivity. To increase the 
number of reference stars, we measured all star images with 
R = 14-21. At this faint end, the star field is very crowded. 
This forced us to improve our procedure of photocentre com- 
putation, which initially was developed for isolated stellar im- 
ages (Lazorenko 2006). We scanned images obtained under var- 
ious seeing conditions and made a detailed census of star po- 
sitions at subpixel precision, of fluxes, and of the point-spread 
function variation across the CCD. This information was used 
to accurately model and subsequently subtract the contamina- 
tion of background counts caused by distant stars. Wings of star 
profiles at distances up to 50 px were approximated by expo- 
nential function with free parameters smoothly varying over the 
CCD. The photocentres x and y were computed by fitting star 
profiles in a 11 x 11 px window with a model with 12 free 
parameters and an auxiliary oscillating function in the central 
5 X 5 px window (Lazorenko 2006). Computations based on the 
Levenberg-Marquardt numerical algorithm of the least squares 
fit were found to give sufficiently good results for stars of ap- 
proximately equal brightnesses at separations larger than about 
10 px. 

3.2. Astrometric model 

Reductions were based on the method previously applied 
to FORS1/2 observations (Lazorenko et al. 2007, 2009). The 
method was shown to efficiently mitigate atmospheric image 
motion, geometric field distortion, DCR, and other effects, thus 
ensuring precision stable over time scales of a few days to a few 
years. For every star in the field, i.e. VB 10 and the reference 
stars, the measured photocenter positions x m and y m in frame m 
at time t were represented by the model 

xq + ^l(x, y)+fi x t + np v + p tan z sin y + d tan z L sin y m 
= x m - ¥(»). k ' 

The left side contains the free model parameters and the model 
function $> hn {x,y). The expression fory-data is similar but con- 
tains cos y instead of sin y. Here, xo is a zero point, k is the mode 
(arbitrary even integer, usually from 4 to 16), <I>jrJ is a polyno- 
mial in x and y of order k/2 - 1 that models the sum of atmo- 
spheric image motion and geometric distortion for each frame 
m. The parameter p x is the proper motion, n is the parallax, and 
p x is the parallax factor in x. The displacement of the star im- 
age due to DCR is modelled by a term with leading parameter 
p, which depends on the star's colour, and containing the zenith 
distance z and the angle y between a direction to zenith and y- 
axis. The next term describes an image displacement opposite 
to that of DCR and introduced by the longitudinal atmospheric 
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Fig. 1. Circular reference fields increasing in size with increasing 
modal number k in the field of view of FORS2 in the case of: a) 
the full data set when an area unavailable on 2 Aug frames was 
cut away at all nights; b) the same but the problem area cut only 
at 2 Aug frames; c) the short data set not using 2 Aug frames. 



dispersion compensator (LADC) of the VLT (Avila et al. 1997). 
This displacement also depends on the star colour via the param- 
eter d « -p. Both p and d are free model parameters and their 
computation does not require external colour data. The LADC 
is automatically adjusted to the average zenith distance zl over 
a given series of frames by setting a distance b between its two 
prisms to b ~ tanzL- Finally, *P(f) represents the induced or- 
bital motion of VB 10 if the planet exists. For reference stars, 
*P(f) = 0. 

Equation 1 defines a system of equations in the combined 
{x, y, rj-domain for which the unique solution is derived under 
the condition that the model parameters of reference stars are 
orthogonal to each basic function of O^. It is solved iteratively 
for all reference stars available in a circular region of radius 
which increases with k and is centered on the target VB 10. The 
optimal field size R^ is the size at which the noise from a refer- 
ence field is equal to the noise from atmospheric image motion. 
The number of reference stars depending on k varies from a few 
dozens to 500. Due to the large number of reference stars, the 
function <t> accurately reproduces the coordinate grid distortion 
introduced by image motion. Solutions for Oa„, at each mode k 
are then used to form equations for the target in the time domain 
only. The solution of this new set of equations yields the tar- 
get's model parameters and position deviations from the model. 
The final output is the average obtained from the solutions at all 
modes k. 

3.3. Comparison of STEPS and FORS2 reference frames 

Because of different zero-points, STEPS proper motion cannot 
be directly applied to FORS2 astrometry. The STEPS reference 
frame is given by 15 bright stars. For the FORS2 reference 
frame, we use both these (except STEPS star Nr. 10 which is a 
binary) and a number of much fainter stars. The system of proper 
motions is defined under the condition £ M;W; = where sum- 
mation is taken over all reference stars ; with weights Wj approx- 
imately proportional to their brightness (Lazorenko et al. 2009). 
In particular, this condition is valid for the proper motions of the 
STEPS stars which, due to their high brigthness (they contribute 
to over 40% of the total light flux), are the basis of the FORS2 
reference frame. Hence, one may expect that Yi^i w il Z w i (the 
weighted mean proper motion in the system of FORS2 proper 
motions) taken over these stars, is nearly zero. The computed 
value of this sum (which is the difference of proper motion zero- 
points) was found to be small: -0.1 ±1.0 and +0.3 + 1.0 mas yr~' 
in RA and Decl. Due to the low precision, we did not apply it but 
included it in the error budget assuming that the uncertainty in 
the proper motion zero point is ~1 mas yr . Thus, both proper 
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motion systems are well consistent, however at time scales of 
~ 0.1 yr only. 

3.4. Short and full dataset 

The data reduction of the first epoch (2 Aug) was problematic 
because of incorrect telescope pointing, owing to which all stars 
in a 300 px wide area just below VB 10 (Fig. la) were imaged 
to another CCD chip and could not serve as reference at this 
particular night. We dealt with this problem in two ways. In the 
first version, we put aside all stars within this area from frames 
of all nights, using only stars outside of this area as reference 
objects. The sizes and shapes of the reference fields in this case 
are shown in Fig. la by five circular segments, each correspond- 
ing to k increasing from six (small radius) to 14 (large radius). 
Alternatively, we used circular reference areas, which for the 2 
Aug night were of the same size but with the lower half vi- 
gnetted (Fig. lb). We also examined a short dataset without the 
2 Aug epoch with circular reference area (Fig. lc). The asym- 
metry of the first two cases degraded the obtained precision in 
comparison to the symmetric configurations of the short dataset. 

We have obtained a solution for all three cases presented in 
Fig. 1. For the final result, we present two solutions. One is ob- 
tained from the short dataset (case V, MJD between 55082 and 
55099 days), and the second one is obtained as the average of 
the results in case 'a' and 'b'. The latter solution is referred to as 
obtained from the full dataset (MJD between 55045 and 55099 
days). 

3.5. Recovering LADC positions 

VB 10 is very red and differs much in colour from the reference 
stars. Therefore, the DCR displacement of this star is very large 
and about 30 mas in Decl and 10 mas in RA. This is taken into 
account by the free model parameters p and d, but their adjust- 
ment requires knowledge of the LADC separation b. 

Information on b was not part of the obtained fits headers and 
was not accessible to us. Therefore we had to solve the inverse 
problem of recovering LADC separations from the observations. 
Because b is fixed for a series of frames of a single night, we 
had to determine a small number of values for b (one value per 
night). Our approach is based on the observation that incorrect 
values of b bias the average night position residuals (x) and (y) of 
field stars. According to Eq. 1, this bias linearly depends on the 
star colour parameter d during one given night. This is illustrated 
by Fig. 2 which shows the distribution and linear dependencies 
of (y) on d for field stars in the case of the short dataset (the 
effect is largest in y). 

Because all observations were made at small hour angles 
within +0.7 h, we initially assumed that LADC was always set to 
the separation corresponding to the meridian. By applying small 
corrections to these initial values of b, we iteratively reached a 
solution without dependence of (y) on colours (lower panel of 
Fig. 2). For these computations, we used all stars within the en- 
tire field of view of FORS2, processing them as target objects 
relative to their own subsets of reference stars. Being based on 
field stars only, this procedure produced a much smaller disper- 
sion of (y) values for VB 10 and for reference stars with large 
colours. Corrections to b were small and within +4% of the 
initial values and were found with precision corresponding to 
+0.025 mas error in (y) for VB 10. For the full dataset, precision 
degraded to ±0.068 mas because we could not use stars in the 
problematic 300 px area below VB 10 (Fig. la). 
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Fig. 2. Distribution of epoch position residuals (y) of field stars 
(dots of different type) as a function of d for four September 
nights and linear dependences of (y) on d (lines of different type) 
for each night. Upper panel: with initial LADC separations b. 
Lower panel: the same with recovered b values. Size of dots is 
proportional to the star brightness. 



3.6. Treatment of colour-dependent terms 

An essential drawback of the above procedure is that the re- 
covery of b introduces small colour-dependent terms p'(d) and 
n'{d), which are similar to proper motion and parallax. This is 
because the above iterations do not garantee convergence to the 
actual values of b ~ tanzL- For example, the restored value of 
tanzL may differ from its actual value by a term which progres- 
sively changes in time as pt tan -1 ziXsiny) , where p is an arbi- 
trary constant and (sin y) is the average value of sin y during the 
given night. Hence the term d tan zl sin y in Eq. 1 generates the 
term dpt sin y (sin y)~ l xdpt linearly dependent on d and which 
therefore can be treated as an extra image motion p'(d) = pd. 
The term p'(d) compensates for the linear change of tanzL in 
time. Similarly, we may assume that the restored values of tanzL 
contain terms proportional to p x and p y . In this case, the solution 
of Eq. 1 for parallax should contain the compensating colour- 
dependent term n'(d). 

When recovering b values, we cannot control the ampli- 
tude of p and of the equivalent parameter ft related to parallax. 
However, they can be detected as proper motion and parallax de- 
pendence on colours (i.e. on d) which we model as a linear trend 
in proper motions of field stars and statistically correct for it. 

The strong correlation between t, p x , and p y does not allow 
us to determine separately p'{d) and n'(d). However, this is not 
required because for short times p x ~ t, we can approximate the 
sum of the proper motion and parallax displacement p'J + n'p x 
by p' x t, where p' is the new effective quantity that substitutes for 
both p' and n'. Thus processing the short dataset, we find and use 
only a single colour term p'(d). The treatment is similar for the y 
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Fig. 3. Orbital motion VP(?) of VB 10 in RA (circles) and Decl 
(triangles) for the full and short FORS2 datasets and the linear 
approximations (solid and dashed lines) of these motions. 



components. Subtraction of p't from the measured positions thus 
simultaneously eliminates the colour term n' (see next section). 
The value of this correction for VB 10 in Decl is p' - 29 ± 
4.5 mas yr~' for the full and p' = 27 + 2.9 mas yr~' for the 
short dataset, respectively. In RA, the corrections are an order of 
magnitude smaller. 

For the full dataset, the above approximation of p.' x t + n' p x 
by p\t is still valid for reference stars, most of which have mod- 
erate d and therefore very small colour-induced parallaxes n' . 
But it is not sufficiently precise for VB 10, because of its large 
colour term d. Consequently, n' is not accurately eliminated by 
applying the correction p'(d) as in the case of the short dataset. 
Therefore, we had to treat the term n' of VB 10 as a free model 
parameter. 

3. 7. Subtraction of parallax and proper motions 

Since the FORS2 observations cover a small fraction (<20 %) of 
the orbital period, VP(f) can be approximated as a linear function 
of t: 



Va(t) ~ Qat, 



(2) 



where Q a = ffV a (f)ldt and Q s = ffV^/dt. With the or- 
bital elements given by Pravdo & Shaklan (2009), we obtained 
G ff (PS) = 13.4 and g 5 (PS) = 22.4 mas yr 1 (Fig. 3). 
Astrometric acceleration terms (deviations from the linearity) 
are smaller than 0.1 mas for the full dataset, thus are negligible. 
For the short dataset, they are even smaller than 0.02 mas. 

In addition, *P(f) is approximately a linear function of the 
parallax factors p x (i) and p y (f), which themselves have an ap- 
proximately linear time-dependence. This causes a strong corre- 
lation between parallax, proper motion, and orbital motion, thus 
makes Eq. 1 degenerated. Therefore, we subtracted parallax and 
proper motion from the measured positions of VB 10. Precise 
values of these parameters were found based on the published 
STEPS astrometric measurements which cover a 9 years long pe- 
riod. The best fit of STEPS data yielded fx a cos(6) = -586.8+0.2 
mas yr -1 , p s = -1361.0 + 0.2 mas yr -1 , and n = 168.0 + 1.2 
mas. These values are very close to the estimates given by 
Pravdo & Shaklan (2009) and Anglada-Escude et al. (2010). 
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Fig. 4. Frame-to-frame deviations in Decl of VB 10 from its 
model motion at zero (VP = 0, filled circles) and non-zero plan- 
etary signal (VP + 0, open circles), the epoch average deviations 
(large diamonds), and the single-frame precisions cr\ (asterisks) 
for the cases 'a', 'b', and 'c' of reference field configurations 
shown in Fig. 1 . 



3.8. Position residuals of models with and without planet 

Subtracting the contributions of parallax and proper motion from 
the measured positions of VB 10 resulted in great simplification 
of Eq. 1, which for VB 10 took the form 



xq + ptanzsiny + litany siny + \n'p x ) 
= x m - p x t - np x - fi'j - oj^(*,y 
yo + p tan z cos y + d tan zl cos y + {n'p v 



(3) 



= y m - p y t - np y - p't - $tlix,y) - ^(f) 



with only four free parameters xq, yo, p, and d in the case of 
the short dataset. For the full dataset, it incorporates the extra 
parameter n' (the term in curly braces, see Sect. 3.6). The right 
side of Eq. 3 contains the measured coordinates x m and y m , the 
in-plate solutions 0^|(x,y) and ®yiXx,y) derived from the refer- 
ence stars, and all other corrections. 

There is a correlation between the orbital signal and the night 
average angles (sin y) which are not zero and tend to increase in 
time. Because of this correlation, the model Eq. 1 filters out any 
component of the signal (e.g. VP), which is linearly dependent on 
t and enters the right side of Eq. 1. Therefore, the output posi- 
tion residuals contain only a part of the initial signal amplitude. 
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Fig. 5. Measured in-frame positions x - <t>' A ', y - <1>'- V ' (squares) 
of VB 10, the model motion (solid line) over the CCD between 
8 and 25 Sept (short dataset) defined primarily by proper motion 
and parallax, and the magnified DCR signature for a single night 
(right upper corner). 



However, this does not hamper our statistical analysis, which we 
performed after the subtraction of *P(?) from the observed po- 
sitions, thus assuming zero input signal and consequently zero 
output signal. In this way, the impact of correlations is mini- 
mized. If the planet does not exist but the subtraction of *P(f) 
was applied, we should detect the inverse signal -¥(?) reduced 
in amplitude because of the correlation between (sin y) and ^{t). 
This is the case corresponding to the last row of Table 2 (Sect. 4), 
where the measured signal (expressed by the parameter Q) has 
about half its expected value -£>(PS). 

Fig. 4 shows the results in terms of model deviations in Decl, 
where the expected signal is largest. The single-frame precision 
0"i includes errors of photocentre measurements, the reference 
frame noise, and the atmospheric noise. It varies from 0.4 to 
0.7 mas depending on seeing. The effect of the vignetted ref- 
erence field of the 2 Aug epoch (configurations 'a' and 'b') is 
seen as a degradation of cr\ to over 1 mas. At other epochs, cr\ 
is larger compared to configuration 'c' because of a larger R^. 

Clear conclusions can be drawn from the short dataset when 
the model (Eq. 3) is most simple and does not require incorpo- 
ration of parallax (Sect. 3.6). We considered the cases with the 
predicted orbital motion subtracted (\P + 0) and with ¥ = 0. 
The epoch average deviations (x) and (y) are very small and ran- 
domly scattered when assuming *T = 0, but display a negative 
trend in time if ^ + 0. Note that small position deviations do not 
correspond to a 'zero' measurement. Instead, they demonstrate 
very precise position measurements, which track the proper mo- 
tion and parallax displacement at the daily rate of 2.7 mas and 
5.0 mas in RA and Decl, respectively. The motion of VB 10 over 
the CCD surface (Fig. 5) for the measurement timespans of 17 
days is 46 mas in RA and 86 mas in Decl, and is dominated by 
parallax and proper motion. DCR effects induce a small-scale 
scatter in the measured positions of one night with an amplitude 
of about 2 mas. Their structure for a typical night is shown in 
5-fold magnification in Fig. 5. 

Similar computations were performed for the full dataset, 
where we had to account for the parallax correction n'(d) (see 
Sect. 3.6), thus used the five parameters xq, yo, p, d, and n' for 
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Fig. 6. Astrometric residuals in the scenario with planet (open 
blue symbols) and without planet (filled black symbols) for the 
short dataset {top panel) and the full dataset (bottom panel) as a 
function of time. Circles and triangles mark residuals in RA and 
Decl, respectively. For clarity, only one errorbar corresponding 
to 0"sum is shown at each epoch, but the uncertainties for the four 
respective measurements are identical. 



fitting of the model Eq. 3. This introduced a strong correlation 
between p x , p y , and *¥, which reduced the amplitude of the de- 
tected signal. 



Table 1. Epoch deviations (x) and (y) in RA and Decl of VB 10 
for the models with Y = and <F + (Eq. 3). 
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T^O 




MJD 

55000+ 


RA Decl 

(mas) (mas) 


RA Decl 

(mas) (mas) 


°"N O" sum CT fl , 

(mas) (mas) (mas) 




Full dataset, 2 Aug-25 Sept 2009 


45.163 
82.032 
92.993 
97.014 
99.011 


-0.26 -0.16 
0.20 0.19 
-0.09 -0.06 
-0.00 0.18 
-0.01 -0.25 


-0.79 0.41 
0.15 0.35 

-0.06 0.11 
0.03 0.08 
0.03 -0.38 


0.19 0.64 0.32 
0.10 0.18 0.16 
0.08 0.11 0.10 
0.07 0.11 0.09 
0.07 0.12 0.10 




Short dataset, 8-25 Sept 2009 


82.032 
92.993 
97.014 
99.011 


0.07 -0.09 
-0.02 -0.04 

0.03 0.07 
-0.06 -0.01 


0.33 0.43 
-0.12 -0.09 

0.02 -0.00 
-0.09 -0.15 


0.11 0.19 0.16 
0.09 0.09 0.06 
0.07 0.08 0.05 
0.08 0.09 0.07 



Table 1 and Fig. 6 summarise the results for the epoch residu- 
als (x) and (y). The astrometric precision is described by a nom- 
inal precision <x N based on errors in photocentre determination, 
the reference frame noise and atmospheric noise. cr sum also in- 
cludes error components which dominate at long time spans and 
originate from the uncertainties in b, p! ', pixel scale, and proper 
motion of VB 10. For the full dataset, it also includes the uncer- 
tainty in the parallax of VB 10. cr^ is the mathematical expec- 
tation of the root-mean-square of (x) and (y), derived from the 
least squares fit (Eq. 3). 
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4. Discussion 

Because the signal was subtracted in Eq. 3, the epoch de- 
viations (x) and (y) should have an expectation value of zero 
and a random scatter irrespective of the model (with or without 
planet). A small dispersion of (x) and (y) values indicates that 
the tested hypothesis is correct. The data shown in Fig. 6 and 
Table 1 is more consistent with *F = (the hypothesis of no 
planetary companion), while incorporating a non-zero planetary 
signal *F + increases the scatter of data points, which does not 
support the existence of the planet VB 10b. We estimated the 
probability P that the planet exists in two ways. 

The first estimate P(x 2 ) was based on the values of x 2 = 
Y,l( x ) 2 / °~lum + (y) 2 /°"sum]' where the sum is taken over epochs. 
The number of degrees of freedom (D.o.F.) is equal to the num- 
ber of epochs minus the number of model parameters (there are 
two and three free parameters for the short and the full dataset, 
respectively). The probabilities P(x 2 ) corresponding to the x 2 
values are given in Table 2. For both the short and the full 
dataset, they are at the equally small at the 2% level if the model 
assumes *T + 0. Because the probability of this event is low, the 
large scatter of (x) and (y) is unlikely to be caused by random 
noise in the observations. 



Table 2. Probabilities P(x 2 ), P(Q) and linear trends Q (mas yr -1 ) 



Dataset 




x 1 


D.o.F. 


P(x l ) 


Q a 


Qs 


P(Q) 


full 





10.6 


7 


0.16 


0.67 


-0.35 


0.36 


full 


*0 


18.0 


7 


0.02 


4.14 


-5.38 


0.02 


short 





1.9 


6 


0.93 


-2.06 


2.47 


0.21 


short 


*0 


14.8 


6 


0.02 


-7.82 


-11.32 


5 ■ 10- 4 



Secondly, we noticed that the observed residuals (x) and (y) 
show a linear trend instead of a random distribution. This can 
be caused by random errors in the observations, but also may 
indicate a wrong value of the signal ¥(?) subtracted in Eq. 3. In 
the latter case, we would expect that the residuals (x) and (y) 
show a linear dependence —Qt opposite to *F(f) (see Eq. 2). 

A large negative trend is seen for the short dataset residuals 
computed with *P + (Fig. 6), whereas they should be near zero 
to support the planetary hypothesis. In contrast, a much smaller 
trend is observed when assuming *P = (no planet). For in- 
stance, the measured values of Q given in Table 2 for the short 
dataset and *P + are -7.82 and -11.32 mas yr~' in RA and 
Decl, respectively. If the null hypothesis \P + is correct, their 
expected values are and mas yr _1 , and -13.4 and -22.4 mas 
yr if the hypothesis is wrong Q¥ = 0). In the last case, the mea- 
sured Q values are always smaller than the values Q(PS) induced 
by the orbital motion, because their magnitude is damped by cor- 
relations between *P(r) and model parameters (cf. Sect. 3.8). 

Because the x 2 criterion is not the most efficient one for the 
characterisation of linear signals, we developed another, more 
powerful approach to obtain an alternative estimate of P. We 
considered Q as quantities describing the dispersion of (x) and 
(y) better than^f 2 and performed Monte Carlo simulations to see 
if the observed features can be explained by random noise in 
the observations. We simulated each data frame by a random 
Gaussian noise with root-mean-square of <r\ and added compo- 
nents modelling the errors in b, and in the proper motion and 
parallax of VB 10. In addition, we included a 1 mas yr~' error to 
account for the uncertainties in the zero point of proper motion 
(cf. Sect. 3.3). 



After fitting the data with the model Eq. 3, we searched 
the residuals for linear trends in time and estimated their co- 
efficients Q* . The false-alarm probabilities P(Q) that a random 
noise produces Q* values larger than the observed ones are given 
in Table 2. In the case of the short dataset and *F + 0, we find that 
the observed linear trends can be explained by Gaussian noise 
with a probability of P(Q) = 5 ■ 10~ 4 , which does not support the 
existence of the planet. For the full dataset, we find P(Q) = 0.02, 
which is not sufficient to draw definite conclusions. In spite of 
the longer timespan and therefore seemingly better conditions 
to characterise a planetary signal, the full dataset does not pro- 
vide a better contraint because of the uncertainty in //, whose 
contribution increases with time. Besides, the planetary signal is 
substantially damped because of its strong correlation with the 
colour correction n' used in Eq. (3) as a free model parameter. 

5. Conclusions 

We conclude that the presence of the announced planet around 
VB 10 is not supported by astrometry. Even assuming half the 
planetary mass (i.e. 3.2 Mj), simulations give a low false alarm 
probability P(Q) = 0.023, which does still not support the ex- 
istence of VB 10b. Our result obtained from astrometry alone 
is in agreement with the conclusion of Bean et al. (2010) and 
Anglada-Escude et al. (2010) based on RV data. 

This study is the first application of the FORS2 camera for 
the search of exoplanets by means of optical astrometry. Because 
of the high astrometric precision of FORS2, the availability of 
external STEPS-based proper motion and parallax of VB 10, 
and the large expected orbital signal, it was possible to perform 
the verification of the planetary companion hypothesis within 
the extremely short observation period of 17 days, which is un- 
usual for astrometric works of this type. The successful use of 
the STEPS data for the reduction of FORS2 observations is justi- 
fied only because we verified that the STEPS and FORS2 proper 
motion reference frames are consistent within 1 mas yr~' uncer- 
tainty. 

We have demonstrated a mean nominal precision of 0.09 
mas per epoch of FORS2/VLT observation for data of reason- 
able quality, despite problems caused by the uncertainty in the 
LADC position. This precision is sufficient for astrometric de- 
tection of planets around ultracool dwarfs. 
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